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We present a perturbative method to estimate the spectral gap for adiabatic quantum optimiza- 
tion, based on the structure of the energy levels in the problem Hamiltonian. We show that for 
problems that have exponentially large number of local minima close to the global minimum, the 
gap becomes exponentially small making the computation time exponentially long. The quantum 
advantage of adiabatic quantum computation may then be accessed only via the local adiabatic 
evolution, which requires phase coherence throughout the evolution and knowledge of the spectrum. 
Such problems, therefore, are not suitable for adiabatic quantum computation. 



It is "widely believed that quantum mechanics can pro- 
vide speedup for certain computations. Different quan- 
tum algorithms have been proposed that potentially can 
solve problems such as factorization llj , unstructured 
search Q , or molecular simulations Q on a quantum 
computer. One type of problems for which quantum me- 
chanics may provide an advantage over classical compu- 
tation is optimization. In optimization problems, one is 
interested in finding solutions that optimize some func- 
tion subject to some constraints. Usually, not only the 
best solution, but also solutions close to it are of interest. 

Physical systems at low temperatures naturally relax 
to their lowest energy states, effectively providing opti- 
mal solutions to their energy function. Such a relaxation 
process, however, may take a very long time. The time to 
reach the low energy states may be significantly reduced 
via an annealing process in which the temperature is re- 
duced from a large value to a small value so slowly that 
the system stays effectively in equilibrium at all times. 
The slow evolution from a thermally disordered to or- 
dered state with decreasing T will settle the system in 
one of its low lying energy states depending on the evolu- 
tion time. Similar ideas have been employed in simulated 
annealing algorithms. 

Quantum annealing (QA) 0, 0] is the quantum ana- 
log of the above classical annealing. In QA the disorder 
is introduced quantum mechanically via a Hamiltonian 
that does not commute with the optimization Hamilto- 
nian. The added term generally has a ground state that 
is a superposition of all the eigenstates of the optimiza- 
tion Hamiltonian. Therefore, the disordered state is a 
superposition rather than a thermal mixture as it is in 
classical annealing. The disorder is removed by slowly 
removing the added term to the Hamiltonian. The sys- 
tem will then settle into one of its low lying energy states 
if the evolution time is long enough. 

Closely related to QA is adiabatic quantum compu- 
tation (AQC)0. In AQC an initial Hamiltonian TLb is 
slowly deformed into a final (problem) Hamiltonian Tip: 



H= [1- s(t)}H B +s(t)H P , 



(1) 



with s(t) changing from to 1 between the initial (£j=0) 
and final (tf) times. In this case, Hb plays the role of 



the disordering Hamiltonian. The main difference be- 
tween QA and AQC is that in the latter, the system is 
constrained to its ground state at all times, starting from 
the ground state of Hb , into which it is designed to be 
initialized, and ending in the ground state of Hp, which 
encodes the solution to the problem of interest. In other 
words, AQC is an exact algorithm while QA is heuristic. 

Unlike QA, AQC is not restricted to optimization prob- 
lems, i.e., the problem Hamiltonian can be non-diagonal. 
For example, a universal AQC can run any quantum algo- 
rithm, and has been shown to be computationally equiv- 
alent to the gate model of quantum computation, as both 
can be efficiently mapped into each other 0, 0] • 

The performance of AQC is determined by the mini- 
mum gap g m between the first two energy levels. In the 
global adiabatic evolution scheme, s is changed uniformly 
with time (s =const.) and the computation time depends 
on g m as r g i b a i oc g~ 2 . In the local adiabatic scheme 0], 
on the other hand, s is a nonlinear function of time de- 
signed in such a way to optimize the computation time 
by spending the majority of the evolution time in the 
vicinity of the anticrossing. As a result, the computation 
time of the local AQC is reduced to ri oca i cx <?~ , which 
scales better with g m . 

The global and local schemes of AQC are also differ- 
ent in terms of their response to decoherence. The global 
scheme is robust against environmental noise and deco- 
herence 13, II , HilTUi, [3] ■ The local adiabatic scheme, 
on the other hand, is very sensitive to decoherence. It 
was shown in Ref. [14| that in order for the local scheme 
to change the scaling of the computation time from oc <?~ 2 
to oc g" 1 , the computation time should be smaller than 
the global dephasing time. Moreover, local adiabatic evo- 
lution requires knowledge of the spectrum which is not 
feasible for general Hamiltonians. 

An important question now is what kind of problems 
can benefit from the quantum advantage of AQC without 
requiring local adiabatic evolution and therefore phase 
coherence? It is known that for the unstructured search 
problem r g i bai = O(N), which is the complexity of 
classical search, while ri oca i = 0(y/~N), which is the opti- 
mal performance of a quantum algorithm. (Here N = 2 n , 
where n is the number of qubits.) Thus, the advantage 
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over classical computation is only possible via the local 
adiabatic evolution. On the other hand, the universal 
AQC 0,101 can provide solution to a problem in poly- 
nomial time if the same problem can be solved in polyno- 
mial time using gate model quantum computation. Evi- 
dently, the polynomial advantage does not depend on lo- 
cal evolution, which can only provide quadratic enhance- 
ment. There have also been previous works to determine 
the complexity of AQC for some other special Hamiltoni- 
ans (15 firL [l8l| . In this letter, we study this problem for 
a rather general form of adiabatic quantum optimization. 

We consider physically realizable initial and final 
Hamiltonians: 



Hp -- 



i—1 — 1 



(2) 
(3) 



where af' z are the Pauli matrices for the z-th qubit, hi 
and Jij are dimensionless local fields and coupling coeffi- 
cients respectively [typically 0(1)}, and £ is some charac- 
teristic energy scale for Tip. The initial Hamiltonian Tig 
has a nondegenerate ground state |-0g(O)) = |0"). Here, 
we have adopted the notation \z) — H® n \z), z 6 {0, l} n , 
for states that are diagonal in the Hadamard basis, with 
H being the Hadamard transformation. 

We denote the ground state of the total Hamiltonian by 
\iPg) — ~^2 z a z\ z )i where a z are complex probability am- 
plitudes. At the beginning of the evolution, a z = 1/y/N, 
therefore \iPg) is a uniform superposition of all the states 
in the computation basis, but at the end of the evolution, 
it is only a superposition of the final solutions. The tran- 
sition from large to small superpositions happens very 
suddenly at the minimum gap, which in the limit g m — > 
represents a quantum phase transition. Here, we only fo- 
cus on first-order phase transition in which the gap is in 
the form of an avoided crossing . If g m is much smaller 
than the separation of the two crossing levels from other 
energy levels, then the slow evolution of the system close 
to the anticrossing will be restricted only to those levels. 
Using a new coordinate e = 2E(s— s*), where E is an 
energy scale characterizing the anticrossing and s* is its 
position, one can write a two-state Hamiltonian: 



H = -{er z + g m T x )/2, 



(4) 



with t x ,z being the Pauli matrices in the two-state sub- 
space. 

Immediately before and after the anticrossing, we write 

fo&±> = hM±co)> = £ af\z), (5) 
ze{o,i}" 

with eo <C E\2, and E%2 being the energy separation 
between the first two excited states. Using (Q} it is easy 



to show that for eo 3> g r 

9m ~ 



eo|(^G>l, 



(6) 



i.e., g m is proportional to the overlap of the wave- 
functions before and after the anticrossing. 
Let us introduce two sets 



{z : \af \ > 6}, 



(J) 



where S is a small number. Since all elements in con- 
tribute to the superposition, the normalization condition 
requires \af\ = 0(l/\/\S ± \), yielding 



N 1 

51 ~7fc 



(8) 



where denotes the cardinality of set S. The minimum 
gap will therefore be 



_\s+ns_ 
7WW 




(9) 



The equality happens when all states in S + also belong 
to S~ . We shall only focus on this case as it provides 
an upper limit for g m . In order to understand what can 
make the gap small, we need to understand how are 
constructed. To this end, we use perturbation expansion. 
We introduce the dimensionless parameter 



s(t)£ 



[!-*(*)] A' 



(10) 



which varies from to oo during the evolution. We shall 
drop the time dependence of £ for simplicity. We start 
by calculating \iPg) near the end of the evolution, where 
( is large, by considering "Hiu = "Hp + (1/C)^s- Let us 
for now assume that the problem has a unique solution, 
therefore the ground state of Hp is a nondegenerate state 

|/). To the 0-th order in 1/C, l^l?) = I/)- Since n B is 
a linear function of erf, it can only generate single qubit 
flips. Thus, to include a state \z) with Hamming distance 
to = \\z — /|| from the solution |/), in \iPg)> we need to 
apply TLb at least to times. Therefore, a f is nonzero only 
after the m-th order perturbation: af=0(l/( m ). This 
restricts the states in S + to be close to / in Hamming 
distance. Requiring l/C" 1 ^^! we find 



S+ w {z: \\z- 
log(l/«5) 



TO C OC 



logC^ 



< m c }, (11) 
C ± ^C(e^±e ). (12) 



The above argument can be easily generalized to 
multi-solution problems by writing the unperturbed 
ground state near the end of the evolution as 



l^ 0) ) 



-1/2 



Yli^i I/')' where fi is the l-th solution 
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among the total N s solutions of the problem. In this 
case, 



S + ps {z : min\\z- fi\\ <m c }. 



(13) 



Therefore, the set S + is constructed from states that are 
close in Hamming distance to the final solutions. 

To find S~ , we perform perturbation expansion around 
Hb, with ( as the small parameter, using Hq = Hb + 
(Hp. Before performing the perturbation expansion, let 
us use our intuition to understand how S~ can be formed. 
At ( = 0, the ground state of the Hamiltonian Hb is the 
uniform superposition of all the states in the computation 
basis. Adding a small perturbation (Hp to the Hamilto- 
nian will introduce a penalty to those eigenstates of Hp 
that have large eigenvalues. As a result, one expects that 
adding (Hp will remove those high energy eigenstates 
from the superposition. The larger the (, the more high 
energy levels will be removed from the superposition and 
eventually only low lying states will survive. 

Let us make this intuitive argument more quantitative. 
Since TIb is diagonal in the Hadamard basis, we need to 
do the perturbation expansion in that basis. Let us write 



(14) 



(0)\ 



ze{o,i}" 

To the 0-th order, the wave function is 
The Hamiltonian Tip is a bilinear function of erf , hence it 
can generate single and two qubit flips in the Hadamard 
basis. Again, in order to include a state \z) into the 
superposition \iPg), where z has a Hamming weight w, 
we need to apply TLb at least w/2 times. This requires 
b z = 0{( w ' 2 ) thereby making b g non-negligible only if 



w < W c (X 



log(VC-)' 



(15) 



In order to determine S~, we need to know how \4>g) 
is formed in the computation basis, not in the Hadamard 
basis. Since Hb does not commute with Tip, one can use 
the uncertainty principle to find the restriction imposed 
by the perturbation in the Hadamard basis, on the wave 
function in the computation basis. Let SEp and SEp 
represent uncertainties in Tip and Hb, respectively: 



SE B = {{Hi) - {Hb?) 1 ' 2 , 
SE P ee {{Hp) - {Hp) 2 ) 1 ' 2 , 



(16) 



where (...) represents expectation value. The uncertainty 
principle requires 5Ep- 5 Ep > \{i\Hb,Hp\) . Every state 
\z) is an eigenstate of Hb with eigenvalue wA, where 
w is the Hamming weight of z. For the ground state 
\ipc), therefore, we have SEb ~ w c A. The uncertainty 
principle requires SEp oc l/w c , leading to 



OC 



{z : E z < E c }, 
log(l/C~) 

io g (i/a) • 



(17) 
(18) 



As expected, the set S is made of low energy eigenstates 
of Hp. 

Equations (Ti"2"f and (jT8j) suggest that as £ — * 1, fn c — > 
oo and E c — > 0. Since the perturbation expansion breaks 
down at ( ~ 1, these equations cannot be extended all 
the way to ( ~ 1- In fact, C ~ 1 is exactly where the 
quantum phase transition and therefore the anticrossing 
occurs. However, to calculate g m using ([6]), we need eo ^> 
g m , which ensures that {ipc) are indeed defined far away 
from the phase transition point, where the perturbation 
expansion and therefore (I13|) and (fT7|) still hold. The 
important fact to notice is that the sets are formed 
in completely different ways: S~ is constructed by all the 
energy levels below some threshold, while S + is formed 
by all the states in Hamming proximity to the answers. 
The two sets could be very different leading to a very 
small energy gap. 

For the upper limit in Q, i.e., g m oc ^|5 + |/|S _ |, 
the computation time will be Ti oca i oc y / Tgi^bIT oc 
^/|iS , ~|/|/S'+|. If the problem Hamiltonian happens to 
have an exponentially large number of low energy states 
that have large Hamming distances to the correct solu- 
tions (i.e., low energy local minima), those states will 
belong to S~ and not to S + . The resulting gap will 
therefore be exponentially small, and the computation 
time will be extremely large. Especially, if \S~\ becomes 
a fraction of A, then r g i b a i = 0(A), which is the com- 
plexity of the exhaustive search. The quantum advan- 
tage then will only be achievable via the local adiabatic 
scheme for which ri oca i = 0{^J~N). 

An interesting example of such difficult problems is 
random 3-Satisfiability problem (3-SAT). These prob- 
lems exhibit a phase transition when the ratio of the 
number of clauses m to the number of variables n reaches 
ss 4.2 20]. Before the phase transition the number of so- 
lutions that satisfy the 3-SAT formula is extremely large, 
but suddenly after the phase transition point the number 
of satisfying solutions drops to zero. Therefore, near the 
phase transition point, by adding a few clauses to the 
3-SAT formula, a large number of states that did satisfy 
it before will no longer do so. Those solutions, how- 
ever, violate only those few clauses. In terms of energy, 
by adding a few terms in the Hamiltonian that provide 
penalties for those few clauses, an exponentially large 
number of states that used to be global minima suddenly 
become local minima but with energies very close to the 
ground state energy. This, therefore, would result in an 
exponentially small gap, as confirmed numerically [21 1 
and analytically [22| . 

Another example is spin glasses [23J, in which the lo- 
cal fields hi are small or zero and coupling coefficients 
randomly couple (only) neighboring qubits. In such prob- 
lems, domains can be formed if a large number of physi- 
cally close qubits are strongly coupled to each other [24( . 
The qubits within a domain minimize the coupling terms 
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Jij in Hp. Those terms, however, remain unchanged if 
all of the qubits in the domain flip together. If the energy 
cost of flipping a domain, imposed by the field terms (hi) 
and by the violation of the bounds (Jij) at the domain 
boundary is not so large, then such a domain flipped 
state will form a low energy local minimum. If the size 
of the domain is large, then the Hamming distance be- 
tween the local minimum and the global one will also be 
large. Thus, the local minimum and all the states close 
to it do not belong to S + , while they do belong to S~ . 
A large number of such domains may make g m exponen- 
tially small. It should be mentioned that if hi = 0, then 
the final Hamiltonian will be symmetric under the total 
spin flip operation and the phase transition is likely to 
be second order, invalidating our assumption. 

To conclude, we have used a perturbative approach to 
estimate the gap size for adiabatic quantum optimization 
problems. The gap is found to be inversely proportional 
to the square root of the number of states that have en- 
ergies close to the global minimum. Therefore, problems 
that have a large number of low energy local minima 
tend to have a small gap. In general, problem instances 
in which the interaction terms in the Hamiltonian domi- 
nate the energy eigenvalues (i.e., typical values of Jij are 
much larger than those for hi) are likely to form low en- 
ergy local minima and therefore make the gap small. If 
the number of low energy local minima becomes expo- 
nentially large, then the gap will be exponentially small. 
In such cases, only a local adiabatic evolution scheme can 
provide quantum advantage over classical computation. 
Local AQC, however, requires phase coherence during 
the evolution and knowledge of the energy spectrum 
which limits its practicality. These problems, although 
unsuitable for AQC, could be suitable for heuristic algo- 
rithms (if approximate solutions are acceptable) , because 
the chance of finding a solution within the acceptance 
threshold will be large. Quantum annealing therefore 
may provide good enough solutions in a short time, al- 
though finding the global minimum via AQC can take an 
extremely long time. 

Finally, it should be mentioned that the energy gaps 
considered here are only the avoided crossing type, which 
correspond to first-order quantum phase transitions. If 
the final Hamiltonian possesses a symmetry that imposes 
a spontaneous symmetry breaking at the anticrossing, 
the resulting phase transition may become second order 



or higher orders. It has been stated that Hamiltonians 
with higher order phase transitions can provide better 
(possibly polynomial) scaling with the number of qubits 
[3, HH]. Study of those is beyond the scope of the present 
paper. 
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